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1. Introduction 

The study of the dependence structure between random variables via copula is a classical problem in 
statistics and other applications. The ease of application of copulas has led to their popularity in various 
areas such as finance, insurance, hydrology and medical studies; see for example. Frees and Valdez (1998), 
Genest et al. (2007) and Cui and Sun (2004). This paper examines the mathematical property of copulas by 
focusing on their lower bound. 

Every copula is bounded by Frechet-Hoeffding lower and upper bounds. While Frechet-Hoeffding up¬ 
per bound corresponds to the maximum copula, Frechet-Hoeffding lower bound is generally not a copula. 
Further, the minimum copula does not exist in general in high dimensions greater than 2; see, for example, 
Kotz and Seeger (1992) and Joe (1997). 

In the insurance and finance field, the maximum copula corresponds to the concept called comonotonic¬ 
ity (Dhaene et al., 2002b). In the respect of risk management, comonotonicity is an important concept, 
because it can be used to describe the perfect positive dependence between competing risks. Importantly it 
provides the solution to various optimization (maximization) problems. However, unlike the perfect posi¬ 
tive dependence, mainly due to the absence of the minimum copula, controversy has remained even in the 
definition of negative extreme dependence In spite of these difficulties, the need for the concept of negative 
extreme dependence has remained in insurance and other applications because it may lead to solutions for 
related optimization problems. Many studies have investigated the negative extreme dependence in various 
contexts. Dhaene and Denuit (1999), Cheung and Fo (2014) and Cheung et al. (2015) defined the concept 
of mutual exclusivity which can be regarded as pairwise countermonotonic movements. On the other hand, 
(Wang and Wang, 2011) proposed the concept of complete mixability, which can be used to minimize the 
variance of the sum of random variables with given marginal distributions. Many papers have recently been 
published in this field (Puccetti et al., 2012; Puccetti and Wang, 2014; Wang and Wang, 2014; Bernard et al., 
2014). While the concepts of mutual exclusivity and complete mixability are both useful in various fields 
of optimization problems, since their concepts both depend on the marginal distributions and are problem 
specific, they may not provide the general concept of negative dependence. 

Fee and Ahn (2014b) proposed a set of negative dependence joint distributions, which is named as 
d-countermonotonic copulas (d-CM). The definition of d-CM is known to be the definition of copula only. 
Furthermore, the set of d-CM copulas is minimal in terms of concordance ordering: there is no copula which 
is strictly smaller in concordance ordering than the given d-CM copula except d-CM copulas. Admitting 
the absence of the minimum element in multivariate dimensions d > 3, the set of minimal copulas can be 
important in optimization problems. However, without understanding the further properties of d-CM copu¬ 
las, choosing the proper d-CM copulas for the given optimization problem can be difficult. Furthermore, as 
specified in Puccetti and Wang (2014), d-CM can be too general to be used for the negative extreme depen¬ 
dence For example, any vector (V, V, ■ ■ ■ , V, 1 — V) with V being a uniform[0,l] random variable is d-CM, 
while it is close to a comonotonic random vector except the last element. Hence in this paper, to remove 
such an almost comonotonic case and emphasize the negative extreme dependence concept, we consider 
only a special subset of d-CM copulas, which will be parameterized by the vector Tit e M+, where is a 
d-dimensional positive Euclidean space. Such set of copulas will be named as 7F-countcrmonoton ic copulas 
(nt-CM). Due to the minimality property of the set of TT-CM copulas, which is inherited from d-CM, we 
expect that the set of 7^-CM copulas might be also useful in various optimization problems. 

However, before we discuss the usefulness of 7it-CM copulas in optimization problems, the existence of 
Tit -CM copulas should be first investigated. While the existence of vt -CM copulas with 

Tit = (1, • • • , 1) e M+ 
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is well known in the literature, see, for example, Lee and Ahn (2014b), existence of 77'-CM copulas is not 
guaranteed for general 7F <G Mi. This paper provides the equivalence condition for the existence of 777-CM 
copulas. For the proof and construction of the copula, we use a simple geometrical method to construct the 
copula. A similar result obtained by using an algebraic method can be found in a recent working paper by 
Wang and Wang (2014). 

Since 777-CM is the property of the copula only, the usefulness of 777-CM may be limited to some 
optimization problems which do not depend on marginal distributions. Puccetti and Wang (2014) also note 
the possible limitedness of 77T-CM (hence ~(1 -CM) in solving optimization problems by commenting that any 
dependence concept which does not take into account marginal distributions may fail to solve optimization 
problems which depend on marginal distributions. Variance minimization of the aggregated sum with given 
marginal distributions, which is formally stated in (19) below, is one such example; detailed literature can 
be found in Gaffke and Ruschendorf (1981); Ruschendorf and Uckelmann (2002); Wang and Wang (2011); 
Puccetti and Wang (2014). As can be intuitively expected, and as will be shown in Section 5 below, it can 
be shown that no single copula universally minimizes the variance of the aggregated sum with arbitrarily 
given marginals. However, we will show that using a set of 777-CM copulas rather than a single copula 
can minimize the variance of the aggregated sum for varying marginal distributions when restricted to the 
uniform distribution family. While our result provides a general solution with no restriction on Tit £ M+, 
a partial solution can be observed in Wang and Wang (2014) for some special cases of Tit £ Mi where 
they are mainly interested in so called joint mixability which aims for the constant aggregated sum. More 
detailed results will be provided in Section 5. 

For a financial application of Trt-CM, we provide a new definition of the herd behavior index. Herd 
behaviors describe the comovement of members in a group. Since herd behaviors in the stock markets arc 
observed usually during financial crises (Dhaene et al., 2012; Choi et al., 2013), measuring the herd behavior 
can be important in managing financial risks. Focusing on the fact that the perfect herd behavior can be 
modeled with the comonotonicity, some herd behavior indices that measure the degree of comonotonicity 
via the concept of the (co)variance have been proposed (Dhaene et al., 2012, 2014a; Choi et al., 2013). 
Measuring the herd behavior using such herd behavior indices can be important as it has been shown to be an 
indicator of the market fear. However, while the concept of comonotonicity is free of marginal distribution 
(and hence so is the herd behavior), these herd behavior measures can depend on marginal distributions, 
as will be shown in Example 2 below. Alternatively, we define the new herd behavior index based on a 
weighted average of bivariate Spearman’s rho. This new herd behavior index is not affected by the marginal 
distributions by definition and will be shown to preserve the concordance ordering. We also show that the 
maximum and minimum of the new herd behavior are closely related with comonotonicity and Tl7-CM. 

The rest of this paper is organized as follows. We first summarize the study notations, and briefly explain 
basic copula theory and countermonotonicity theory in Section 2. The concept of 77T-CM is introduced in 
Section 3, and the existence of 77'-CM copula is demonstrated in Section 4. Section 5 applies the concept of 
TT^-CM to variance minimization problems. The definition and minimization of the new herd behavior index 
are discussed in in Section 6, which is followed by the conclusions. 

2, Notations and Preliminary Results 

2.1. Conventions 

Let d > 2 be integers and W 1 denotes d dimensional Euclidean space. Especially, let be d dimen¬ 
sional positive Euclidean space. Further [a, b] x [a, b] ■ ■ ■ x [a, b\ C U. d is denoted by [a, b] d . We use to 
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denote d-variatc vectors: especially, lower case 


x = (x±,x 2 , ■■■ ,x d ) 

denotes constant vectors in W 1 and upper case 

l = (X l ,X 2 ,--- ,X d ) 

denotes d-variatc random vectors. More specifically 

u := (ui, ■ ■ ■ ,u d ) and it := (wi, ■ ■ ■ , w d ) 

will be used to denote constant vectors in [0, l] rf and respectively. Finally, use V to denote a uniform [0, 1] 
random variable. 

Unless specified, we assume A^ be a d-dimensional random vector having H as its cumulative distribu¬ 
tion function defined by 

H(lt) = P(Xi < xi, ■ ■ ■ , X d < x d ) for it £ R d , 

and the marginal distribution of X, is Fi(y) := P(Xi < y) for i G {1, - - - , d} and y 6 R. Define 
J- d (F \, • • • , F d ) to be the Frechet space of d-variatc random vectors with marginal distribution F\, ■ ■ ■ , F d . 
Hence, ^ € F,i(F\. • • • , F d ). Equivalently, we also denote Ft £ Fd(F \, • • • , F d ). We use T d to denote the 
special case of Frechet space, where all marginal distributions are uniform[0,1], 

This paper assumes that marginals distributions arc continuous. According to Sklar (1959), given H £ 
F(Fi, • • • , F d ), there exists a unique function C : [0, l} d —> [0,1] satisfying 

H(lt) = C(Fi(xi), ■ ■ ■ ,F d (x d )). 

The function C is called a copula, which is also a distribution function on [0, l] d . Further information on 
copulas can be found, for example, Cherubini et al. (2004) and Nelsen (2006). 

Any H £ F(Fi, • • • , F d ) satisfies 

WiFxix!),--- ,F d (x d )) <H(lt) - ,F d (x d )), for all it £ M d , 


where 

W(lt) := maxjrti + • • • + u d — (d — 1), 0} and M(~it ) := min{iii, • • ■ , u d }, (1) 

for It £ [0, \} d . W and M in (1) arc called the Frechet-Hoeffding lower and Frechet-Hoeffding upper 
bounds, respectively. Note that M(F \, • • • , F d ) is a cumulative distribution of a d-variatc random vector 
while W(F\. ■ ■ ■ , F d ) is not in general. Let id be a survival distribution function defined as 

Htf) :=P(X i > xi, ■ ■ ■ , X d > x d ) for it £ R d . 

For H , H* £ F(F\, • • • , F d ), the concordance ordering H A H* is defined by 

H(lt) < H* (it) and H(lt) < H*(lt) for all it £ R d . 
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Furthermore, define H = H* if 


H(lt) = H*(!t) 

for any it G M rf . Equivalently, denote It = X* if H = H*, where the cumulative distribution function of 
X* is H*. Unless specified, 

#:=([/!,••• ,U d ), IT :=(Ut,- •• ,u* d ), and U* := (U **, • • ■ , U* d *) 


arc ci-variatc random vectors 
respectively. For example, 


in X,i having copula C, C* and C** as cumulative distributions functions. 


P (Ui < ui ,-- • , U d < u d ) = C(lt) 


for it G [ 0 , l] d . 

It will be convenient to define the minimal and minimum copulas. For d > 2, we define (/-dimensional 
copula C G Td as a minimum(maximal) copula if the inequality 


C* >- {X)C 


for any d-dimensional copula C* G J 7 ,/. Similarly, for d > 2, define d-dimensional copula C G T,i as a 
minimal(maximal) copula if the inequality 

C* < (y)C 

for some d-dimensional copula C* G T,i implies C* = C. Define the set of copulas C C T ( i to be minimal 
in set concordance ordering if any C G C and C* G J- C i with 

C* 


implies 


C* G C. 


By definition, C is minimal in set concordance ordering if C is empty. Clearly, the definition of minimality 
in set concordance ordering is a weaker concept than the definition of minimal copula. In the minimality of 
set concordance ordering, the quality of the minimality depends on the size of the set. For example, Frechet 
space is minimal in set concordance ordering. On the other hand, if C has a single element, the definition of 
the minimality in set concordance ordering coincides with the definition of the minimal copula. 


2.2. Review of d-Countermonotonicity 

Comonotonicity has gained popularity in actuarial science and finance. Conceptually, a random vector 
it is comonotonic if all of its components move in the same direction. Comonotonicity is useful in several 
areas, such as the bound problems of an aggregate sum (Dhaene et ah, 2006; Cheung and Vanduffel, 2013) 
and hedging problems (Cheung et al., 2011). Recently, comonotonicity has been used in describing the 
economic crisis (Dhaene et ah, 2012, 2014b; Choi et al., 2013). 

Countermonotonicity is the opposite concept to comonotonicity. Conceptually, in the bivariate case, a 
random vector it is countermonotonic if two components move in the opposite directions. The following 
classical results summarize the equivalent conditions of countermonotonicity in bivariate dimensions. 

Definition 1 . A set A C M 2 is countermonotonicf comonotonic) if the following inequality holds 

(afi - yi){x2 ~ 2/2) < (>)0 for all it, ft G M 2 . 
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A is called countermonotonic(comonotonic) if it has countermonotonic(comonotonic) support. 
Theorem 1. For a bivariate random vector X, we have the following equivalent statements. 

i. is countermonotonic 

ii. For any 

p(2<-2) = max (F 1 (xi) + F 2 (x 2 ) - 1, 0} 

Hi. For Uniformed, 1) random variable U\, we have 

l±{Fi\Th),Fi\\-lhj). 


( 2 ) 


While the extension of comonotonicity into multivariate dimensions d > 2 is straightforward, there is no 
obvious extension of countermonotonicity into multivariate dimensions d > 2. As discussed in Lee and Ahn 
(2014b), the difficulty of the extension of countermonotonicity arises partially due to the lack of minimum 
copula. In this paper, we provide a set of minimal copulas, which can be viewed as a natural extension of 
countermonotonicity in two dimension into multivariate dimensions. 


3. Weighted Countermonotonicity 

As an extension of countermonotonicity or negative extreme dependence in multivariate dimensions, 
Lee and Ahn (2014b) introduced the concept of d-CM. While d- CM copulas arc theoretically interesting, 
the existence and construction of d-CM copulas with certain parametric functions remain unknown, and it 
may therefore be hard to apply d-CM copulas to various optimization problems. Furthermore, the concept 
of d-CM may be too general to describe the negative dependence concept as briefly specified in Puccetti and 
Wang (2014), where the example of (V.V. - ■ ■ , V. 1 — V) was given. Alternatively, Lee and Ahn (2014b) 
introduced the concept of strict d-CM as a special case of d-CM, which is useful in some minimization 
problems. However, because of the symmetricity of strict d-CM, it cannot be used for non-symmetric opti¬ 
mization problems, as will be explained in Section 5. For completeness in the paper, we have summarized 
the definitions and properties of (strict) d-CM in the Appendix. 

In this section, we introduce a new class of extremal negative dependent copulas, which wifi be called 
77’ -Countermonot(>n ic (ut-CM) copulas and can be interpreted as a set of minimal copulas as shown in 
Corollary 1 below. Remark 1 addresses that the set of A-CM copulas can be interpreted as generalized 
strict d-CM, and further shows that 7^-CM copulas are the subset of d-CM copulas. 

Definition 2. A d-variate random vector X is X-CM if 


p(j2mF i (x i ) = ^f^ 


= l. 


s i=l 


Equivalently, we say that H is vft-CM if 7? is v$-CM. Furthermore, when 7? is -CM, we define ~v$ as a 
shape vector of 7?. 

nt-CM can be regarded as multivariate extension of countermonotonicity into multivariate dimensions. 
First, assume that 2? G F 2 (F\, F 2 ) is countermonotonic. Then, since A is a continuous random vector, 
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Theorem 1. iii concludes that 

Pi(*i) + F 2 (X 2 ) = F 1 o Fi\U!) + F 2 o Ff\ 1 - C/i) 

which in turn implies 3? is T-CM for any wm = w 2 > 0. On the other hand, assume that A' is F-CM with 
wi = w 2 > 0. Then by Definition 2, we have 

F 1 (X 1 ) + F 2 {X 2 ) = 1 

with probability 1, which in turn concludes that the support of A^ is countermonotonic. So we can conclude 
that X £ F 2 (F\, F 2 ) is countermonotonic if and only if A" is F-CM with mi = w 2 . 

As can be expected from Definition 2, is a property of copula only, and this is summarized in the 

following lemma. The proof is si mi lar to that of Lemma 1 in Lee and Ahn (2014b). However, the result in 
the following lemma is more useful as it shows that the shape vector is invariant to marginal distributions. 

Lemma 1. Let X and X* be random vectors from the distribution functions 

H = C(F\, • • • , Ff) and H* = C(FF ■ ■ ■ } F%), 


respectively, where marginal distribution functions, F\, ■ ■ ■ , F r j, are possibly different from marginal distri¬ 
bution functions, Ff, ■ ■ ■ ,Ff. Then X is F-CM if and only if XT is F-CM. 

Proof Since two random vectors (Fi(Xi), ■ ■ ■ , F,i(X c j)) and ((A"(), • • • , Ff ( X* d )) have copula C as 
the same distribution functions, we have 




if and only if 





Hence we conclude that X is ri^-CM if and only if X* is 


= 1 , 


= 1. 


( 3 ) 


( 4 ) 

□ 


In the following definition, we provide the copula version of TT-CM. Note that, for the property of F- 
CM, it is enough to study the copula version of nt-CM, because nt-CM is a property of copula only. Hence, 
throughout this paper, we will use the following definition as the definition of r^-CM. 

Definition 3. A d-variate random vector j} is F-cm if 


pi^WiUi 


E d 

i=i w i 


= l. 


( 5 ) 


\i=1 


Equivalently, we say that C is F-CM ifTf is F-CM. Here, F is called as a shape vector of A. 
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Remark 1. Note that 77-CM is d-CM with parameter functions 

fi(y) = CiFi{y) 


for i £ {1, • • • , d\ and y £ M, where 

2 Wi 

c * := ~d -' 

E wj 

3 = 1 

Furthermore, since Tv-CM coincides with strict d-CM when w\ = • • • = w d , the set of strict d-CM copula 
is the subset ofi7-CM copulas, in the Appendix. For convenience, we summarize the definitions of d-CM 
and strict d-CM in Definition 6 and Definition 7 in the Appendix. 

The following corollary explains that the set of T-CM copulas can be regarded to have minimality in 
concordance ordering as a set. Since nt-CM is a special case of d-CM as shown in Remark 1, the proof of 
the following corollary is immediate from Lee and Ahn (2014b). However, for completeness in the paper, 
we present the proof in the Appendix. 

Corollary 1. For given Tfi £ let C be the set of 7$-CM: i.e. C is defined as 

C := {C £ T d | C is T7-CM) . 

Then C is minimal in set concordance ordering. 

As briefly mentioned in Section 1, since there is no minimum copula available for d > 3, it is clear that 
minimal copulas will play a key role in various minimization problems. In this sense, Corollary 1 addresses 
an important property of n^-CM copulas: the set of T-CM copulas achieves minimality in the sense that 
there are no copulas strictly smaller than the i^-CM copula other than 77T-CM copulas. Hence, the concept 
of -u^-CM can be useful in various minimization/maximization problems as will be explained in Section 5 
below. For a discussion of the usage of ut-CM copulas, it is essential to check the existence of r^-CM 
copulas as will be shown in the next section. 

4. Condition to Achieve the Weighted Countermonotonicity 

Depending on the given marginal distributions, (5) may not be always achieved. For example, for 
(wi, wf) = (2,1), none of (Ui, U- 2 ) £ F 2 can achieve the condition in (5). In this section, we provide the 
equivalence condition of the weight If £ for the existence of FF-CM copula. Note that a similar result 
can be found in a recent working paper by Wang and Wang (2014), which explains an algebraic way of 
constructing nt-CM copulas. We first define the set of weights where FF-CM copulas exists. 

Notation 1. Define the set of weights in 3 -dimensions as follows 

f 3 

W 3 := < (wi,w 2 ,w 3 ) £ M+| y ^Wj > 2max{roi,W2,w 3 } 
l i=i 

Note that the set W 3 is equivalent with the set of the line lengths in triangles (including degenerate 
triangles). 

Lemma 2. For any 7$ £ W 3 , there exists Tf.-CM copulas. 




Proof. For convenience, define 


zi : = 


W2 +W3- W 1 W 3 +W1-W2 , VJl + W2 - w 3 3 

and z 3 := ---, for some (■ w\,w 2 ,w 3 ) G Mi 


22 := 


2 w 3 


2W2 

and denote it E if 

Ul W\ +U2W2 + u 3 w 3 = 

Now, let us consider the following three points 


2iu\ 

W 1 +W 2 + w 3 


( 6 ) 


pt := ( 1 , 21 , 0 ), p 2 : = ( 0 , 1 , Z 2 ) and pt := (z 3 , 0 , 1 ), 

and observe that the points satisfy pf G 'I>( wt j for i G {1, 2,3}. Hence any point on the line that connects 
pi and pj is again in <J>( u{): i.e. 

tpl + (1 - t)pj G (7) 

for any 0 < t < 1 and i,j G {1,2, 3}. Further, by the assumption Tit G W 3 , the following inequalities can 
be derived 

0 < Zi < 1 for i G {1, 2, 3}, 

which in turn implies 

tPi + (1 - t)pj G [0,1] 3 (8) 

for any 0 < t < 1 and i, j G {1,2, 3}. Note that the trace of ( 8 ) is triangular in [0, l ] 3 with vertices lying on 
Pi, P 2 and p|. 

Now for the given triangle with vertices pt, pt and p 3 , we give positive weights mi . 2 , m 2 , 3 , m. 3,1 to 
each edge p\p 2 , P 2 P 3 and p 3 p\ such that weights arc uniformly distributed on each edge. Here we assume 
that that the sum of weights is given as mi ,2 + m 2,3 + m. 3,1 = 1 , so that the weights on the edges of the 

triangle define a random vector A" = (X\. X?, X 3 ). In defining H as the cumulative distribution function 

of the random vector our goal is to show that there exist the weights (mi, 2 , m 2 , 3 , 1713 , 1 ) G M'{_ which 
make H a copula. 

To show that H is a copula, it is enough to show that H is 2-increasing and that the marginals of H arc 
a uniform[0,1] distribution (Nelsen, 2006). Since H is defined by the nonnegative weights m. 12 , m 23 , mm 
that arc distributed on the edges of the triangular, it is obvious that H is 2-increasing. Now, it remains 
to show that marginals of H are a uniform[0,1] distribution. Since weights mi 2 - 77 / 23 , ni 3 \ arc uniformly 
distributed on each edge, it is enough to check uniformity on each vertex of the triangle, which is equivalent 
to show 

P (Xi < z 3 ) = zi, P (X 2 < z 2 ) = z 2 and P (X 3 < z 3 ) = z 3 . (9) 

Each equation in (9) is equivalent with 

Z; 1 = m i2 • 0 + 777-23 Z\ + m 3 l • 1, 

Z 2 = m \2 ■ 1 + 77723 • 0 + m 3 \ z 2 , (10) 

Z 3 = mi 2 z 3 + m-23 • 1 + mm • 0, 
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respectively. The solution of (10) is 


m 12 = 


m23 = 


m 31 = 


(1 - Zl)(l - Z 2 )( 1 - z 3 ) - (1 - zi)(l - 2 : 2 ) + (1 - zi) 
(1 - zi)(l - z 2 )(l - z 3 ) + 1 
(1 - zi){l - Z 2 )( 1 - Z 3 ) - (1 - ^2)(1 - Z 3 ) + (1 - Z 2 ) 
(1 - zi)(l - z 2 )( 1 - Z3) + 1 
(1 - zi)(l - Z 2 )(l - Z 3 ) - (1 - z 3 )(l - Zi) + (1 - z 3 ) 
(1 - Zi)(l - z 2 )( 1 - Z 3 ) + 1 


While + 7772,3 + 7 /( 3 , i f 1 for general (z\. z 2 . Z3) G [0, l] 3 , a tedious but straightforward calculation 
shows that, with (zi, z 2 , Z3) defined in (6), (7711,2, 7772,3,7773,1) defined in (11) always satisfies 


7771,2 + 7772,3 + 771.3,1 = 1. 


Finally, (7) derives that + is 7?-CM. 


While (zi, Z2, Z3) is some vector in [ 0 , l] 3 , it is worth mentioning that the definition (6) is crucial to 
guarantee that the solution (771-1,2,7772,3,7773,1) of ( 9 ) satisfies 

7771,2 + 7772,3 + 771-3,1 = 1. 

In other words, (771-1,2,7772,3,7773,1), which satisfies ( 9 ) may not satisfy 

771-1,2 + 7772,3 + 777 3 ,l = 1 

for and arbitrarily given (zi, z 2 , Z3) <G [ 0 , l ] 3 that does not satisfy the condition ( 6 ). For example, for the 
arbitrarily given (zi, z 2 , Z3) = ( 0 . 5 , 0 . 3 , 0 . 2 ), the solution (7771,2,771-2,3,7773,1) of ( 9 ) defined in ( 11 ) has 

7771,2 + 7772,3 + 771-3,1 > 1. 

The following lemma is an extension of Lemma 2 into multivariate dimensions d > 3 . 

Lemma 3. For given £ M+, if there exist disjoint subsets A, B and C of {w i, • • • , vj,i } such that 


Wi , Wi, yy Wi I £ W 3 and Au B U C = {tui, • • • , Wd] , 

Ui^A Wi^B WidC ] 


then there exists a random vector U whose marginals are uniform[0 , 1] and it satisfies 


y>^ 


E d 

i=i w i 


Proof Let lJ be random vectors with marginals being uniform[0,l]. Further, let 

Ci := £ Ui , C 2 := E Ut and C 3 := J2 Ui ~ 
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Now the proofs are trivial if we set Ufs in the same subset as being comonotonic i.e. U, and Uj arc 
comonotonic if either i, j £ A, i,j £ B or i.j £ C. □ 


The following lemma provides the equivalence condition of (12), which is more intuitive and easy to 
verify. 

Lemma 4. For the given weight ~v5 £ R(j_, we have the following inequality 

d 

max{w L • • • , w d } <^Wi- max{w b • • • , w d } (14) 

i= 1 

if and only if there exist disjoint subsets A, B and C of {w \, • • • , w d } satisfying (12). 

Proof First observe that (12) implies (14) is trivial. Hence it remains to show (14) implies (12). Without 
loss of generality, let w\ > • • • > w d . For any integer d > 3, define 

w *2 := ^2 Wi and w 3 := Wi 
i&iO itz^E 

where Zo := {wi\i f 1 and i < d is odd number} and Z 77 := < d is even number}. Then it is 

straightforward to show that 

vj] +W 2 >u >3 and w\ + w\ > wZ, (15) 

Hence, along with (15), if we assume that ul satisfies (14), we can conclude (w\, w^w^) G W 3 , which in 
turn implies (12) with A = {mi}, B = {wi\i G Z 0 } and C = {wi\i G 1>e\- □ 

So far in Lemma 3 and Lemma 4, we have provided sufficient conditions for the existence 7T-CM 
copula. Then the natural question is to check whether they arc also necessary conditions or not. The 
following corollary shows that the condition in (14) is also a necessary condition for the existence of r^-CM 
copulas. 

Corollary 2. For the given weight IF G M+, there exists random vector IF -CM random vector U if and 
only if dF satisfies 

d 

max{uq, • • • , w d } < ^2 Wi — max{roi, • • • , w d }. (16) 

i —1 

Proof It is enough to show that t^-CM implies (16). First, consider a weight (w \, • • • , w d ) G Ml such that 
one weight, say w \, is greater than the sum of all other weights 

d 

W 1 > ( 17 ) 

i=2 

Then, it is obvious that there does not exist any random vector l2 whose marginals arc uniform[0,1] and 
satisfies (5): this can be easily verified using the following variance comparison; 

/ d 

War (w\U\) > War I ^ WjUj 
\i =2 
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Hence, we can conclude that there does not exist Trt-CM copula under the condition (17), which concludes 
the claim. 

□ 



Figure 1: w-CM Copula with (wi, m 2 , m 3 ) = (5,4, 3) 


Remark 2. For any Tit £ Ml satisfying (16), the choice of Oft-CM copula is not unique. For example, 
in the proof of Lemma 2, we show how to construct Tit-CM copula C for d = 3. On the other hand, the 
construction method in the proof of Lemma 2 and the following choices of z *, z%, z\ £ [0,1] defined as 


* W 2 + m 3 - wi * m 3 + mi — m 2 , * mi + m 2 - m 3 

z 1 := -, := - and z-, := -. 

1 2w 3 ’ 2 2wi 3 2 w 2 

and three points pf, pf, pf £ [0, l] 3 defined as 


pi := ( 1 , 0 , zi), ^ := (z 2 , 1 , 0 ) and p 3 := ( 0 , z 3 , 1 ), 


vv;7Z derive another choice of -CM copula C* with C f C*. 

In the bivariate case. Corollary 2 concludes that (Z7i, Uf) being 77 r -CM implies that 10 \ = W 2 which 
coincides with the concept of countermonotonicity as we already mentioned in Section 3. The following 
example shows the numerical example of the construction of 77-CM copula using the logic in the proof of 
Lemma 2. 


Example 1. Let (mi, m 2 , m 3 ) = (5,4,3). Since 


5 = mi < m 2 + m 3 = 4 + 3, 
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we know that ~v5 G W 3 and, by Corollary 2, there exists a ~v$-CM random vector U G J- 3 . Using the 
techniques used in (2), we can construct -CM random vector having mass m± 2 , m .23 and uniformly 
distributed on the each edge p\p\, P 2 P 3 andpfpi, respectively. Here 


and 


Pi = U’T’ 0 > P 2 = 0,1,- and p 3 = —, 0,1 


6 3 2 

m 12 = —, m 2 3 = — and m 3 1 = — 


Hence, for example, we have 

C (1,0.25,1/3) = 1/3 * m 3 i 
2 

“ 33' 

Finally, Figure 1 shows the support of random vector JJ. 


(18) 


5. Application to the Variance Minimization Problem 

Finding the maximum and the minimum of variance in the aggregated sum with given marginal distri¬ 
butions is the classical optimization problem 

Var ’ for given Xj ~ Fj, i = l,---,d. (19) 

First of all, the maximization of (19) is straightforward using the comonotonic random vectors. For the 
minimization problem with d = 2, the answer is trivial with countermonotonic random variables. Regarding 
general dimensions d > 3, minimization of (19) was solved for some cases of marginal distributions (Gaffke 
and Riischendorf, 1981; Rilschendorf and Uckelmann, 2002; Wang and Wang, 2011; Puccetti and Wang, 
2014). However, minimization of (19) is not easy in general for d > 3. The following remark, which can be 
easily derived from Theorem 2.7 of Dhaene et al. (2014b), states that variance minimization problems are 
related with concordance ordering, which may offer some hints in the minimization of (19). 

Remark 3. Let F \, • • • ,F d be distribution functions having finite variances. If 

X 1 , leF d (F i,--- ,F d ) 


with C* F C, then 


d \ / d 

E - Y .* £ Var Or X, 

i=l J \i 1 

From the remark, it is clear that the minimization and maximization of (19) is related with a minimum 
and maximum copula. While maximization of (19) is related with the comonotonic copula, due to the 
absence of the minimum copula for d > 3, the minimization of (19) is related with the set of the minimal 
copulas. Of course, the choice of the proper set of minimal copulas depends on the marginal distributions. 
Among many other choices of the marginal distributions in (19), this paper considers the uniform marginal 
distributions as shown in the following definition, which may be the simplest versions of (19). The following 
assumption is useful to simplify the notation in several theorems in this section. 
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Assumption 1. Assume that w\ = maxjiyi, • • • , w d }• 
Definition 4. For given Tit G M+, define 




(Tlt) := inf \ Yar ( ^ U t 


vi=l 


Ui is uniform ]}), m*] random variables, i = 1 , • • • , g? > (20) 


and 


m+(Tit) := sup < Var ( (T* 


vi=l 


Ui is uniform [0, Wi] random variables, i = 1, • • ■ ,d 


where Ui is uniform [0, m*] random variables for i = 1, • • • , d. 
Equivalently, m_(Tit) and rn + {~if) can be written as 


m_( Tit) = 


inf | Var G j 


and 


The upper bound 


m_i_(nt) = sup I Var WjU^j 




G -TT 


1 / d 

? n+(TZt) = — ^ 




U=1 


is achieved if and only if lJ is comonotonic (Kaas et al., 2002; Dhaene et al., 2002a,b). Regarding the lower 
bound, when 


2 max{wi, • • • ,w d } < £ Wi, 


( 21 ) 


i=1 


Corollary 2 concludes that 

m_(nt) = 0. 

However for Tit which does not satisfy (21), minimization is not straightforward. 

For Tit G Mtj. which does not satisfy (21), Theorem 2 below finds the explicit expression for m_(Ht). 
More importantly, we also show that the minimum m_(Tlt) is achieved with nr-CM copulas even though 
v? G M([ may not be the same as Tit. Finally, Corollary 3 provides the complete solution for m_(Tlt) 
and m + (Ht) for any given Tit G M([. Before we examine the main results, it is convenient to present the 
following lemma and notations. 

Lemma 5. Let Tit G Ml satisfy Assumption 1 and 


d 


2 U ’1 = Y1 


= > Wi . 
i=i 


( 22 ) 
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Then the following inequality holds 


cov 


Ui,^WjUj 


i— 1 


>o, 


(23) 


where the equality holds if and only iflJ is vS-CM. 
Proof. We first prove the inequality (23) as 


cov 


U\,'y"w l U l 


i =1 


= y^cov [Ui,WjUi\ 
i =1 

d 

= wfVar (Ui) + Y mcovr [U u Ui] y/VariUjy/VariUi) 

i =2 
d 

> wfV ar (Ui) — £ Wiy/V ar (Ui)s/Var (Ui) 


i=2 


d 

w\ - Y w i 

i =2 


Var (Ui) 


= 0 


where the inequality arises from the fact that correlation of any two random variables is greater than — 1, 
and the last equality is from the condition (22). Furthermore, since nt satisfies the condition (16), Corollary 
2 concludes that the inequality in (23) becomes equality if and only if Tf is t^-CM. □ 


Notation 2. For given IF e M'|, define 


W; 


w* := (wt, ■ ■ ■ 


> w d) as 


d 

WT, 

if 2 Wi< Wj 

3= 1 

d 

d 

JfWi- wfi 

if 2 Wi > wj 

J =1 

3 = 1 


for i = 1, • • • , d. Then, one can easily confirm w £ R(j_. Further, let 


(24) 




d \ 

2 max{wi, • • 

•, m} - Y Wi ) 


i=1 / + 


Since there always exists t^-CM random vector U G T,\ for TF £ M'|, which satisfies the condition 
(16), we conclude rri-ifw) = 0 in this case. The following proposition provides the tight bound of m-(lF) 
when iF does not satisfy the condition (16). The main idea of the proof is to shrink the largest weight so 
that the new weights satisfy the condition (16), which in turn results in constant summation or zero variance. 
Then Lemma 5 shows that only the remaining paid of the largest weight contributes the lower bound of the 
variance specified in (20). 
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(25) 


Theorem 2. Let the weight vector IF £ Ml satisfy Assumption 1 and 


2 w i > ^ 


Wi. 


i= 1 


Then the following inequality holds 


l(v$) < \ar WjU^j 


(26) 


where the inequality is attained if and only if l) is w*-CM. Furthermore, for any iF £ satisfying the 
condition (25), there exists random vector j} £ fFd which achieves the equality in (26). 

Proof. First, observe «'i — w) > 0 for the given IF satisfying the condition (25). Further we have 

2w\ = 2 f Wi — 


v i=l 


K 1=1 


i =2 


^2uh-w i ) 
i- 
i 

E 


Wi 


(27) 


= z^ w 

i =1 


Since w\ — w\ > 0, we have 


Var ^ Wi u)j = Var ((mi - wfil 7i) + Var (w*^ + w i U ^j 

d 

(wi-w\)JJi, wfUi + 'y^ j w i U i 


+ cov 


i =2 


Var ((m - <)(7i) + Var ^ m*f/^ 

d 

Ui, wllh + JZ'U’iUi 


(28) 


i=2 


+ (mi — wl) COV 

> Var ((mi — m*)(7i) 

= ^(wt - O 2 , 

where the last inequality is from (27) and Lemma 5. Furthermore, since variance of the random variable is 
0 if and only if the random variable is constant with probability 1, Lemma 5 concludes that the inequality in 
(28) is equality if and only if T? is w*-CM. 

Finally, since w* satisfies the condition (27) (hence the condition (16)), there always exists w*-CM 
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random vector if G J- d , which in turn implies that the equality in (26) can be always achieved. 


□ 


Based on Theorem 2, the following corollary provides the complete solution for the optimization prob¬ 
lem in Definition 4. 

Corollary 3. For the given IF £ M+, the following inequality holds 

l(u$) < \ar • 


v£=l 


(29) 


\i =1 


The lower bound of (29) is attained if and only if if is uf-CM and the upper bound of (29) is attained if 
and only if U is comonotonic. 

Proof. The upper bound in (29) is a classical result which can explained by comonotonic if ; see Kaas et al. 
(2002); Dhaene et al. (2002a,b) for details. Hence it is enough to show the lower bound 
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d 


1 2 


2 max{u>i, • • • ,w d } - ^ 


~ 2^ Wi 
i =1 


+ 


< Var ( ^WjUj 


(30) 


s. 2—1 


and the equality holds if and only if if is w*-CM. 

For the lower bound in (29), consider two cases depending on IF. First, consider the following condition 


on w £ 


2 max{wi, • • • , w d } — Wi < 0. 


i —1 


Since l(lF) = 0 in this case, nonnegativeness of the variance shows the left inequality of (29). Furthermore, 
Corollary 2 implies that the equality in (29) holds if and only if JJ is u/ -CM. Finally, for lF satisfying 

d 

2 max{wi, • • • , w d } - ^ w l > 0, 
i =1 


the same result was summarized in Theorem 2. □ 

6. Marginal Free Herd Behavior Index 

In this section, we define the marginal free measure of dependence which can be interpreted as the 
measure for the herd behavior. We first review various herd behavior indices, which measure the degree of 
comovement or comonotonicity (Dhaene et al., 2012, 2014a; Choi et al., 2013). While such indices need 
to be marginal free (because the concept of comovement or comonotonicity is a definition of copula only), 
in Subsection 6.2, we observe that such indices can be distorted by marginal distributions. Alternatively, 
Subsection 6.3 presents a definition of measures of dependence that is free of marginal distributions. 
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6.1. Review of the Measures of Dependence 

Herd behavior is a general concept often used in various fields such as financial and psychology to 
describe the irrational comovement of members in a group. The recent financial crises have further high¬ 
lighted the importance of understanding the herd behavior. There have been several attempts to measure the 
herd behaviors through herd behavior indices. In this subsection, we will briefly review some known herd 
behavior indices in the financial context. 

Let .V be d individual stock prices at a time t assuming that the current time is fixed at 0. For the given 
X, the market index S is defined as the weighted sum of the d individual stock prices: 

d 

s = wjXj, 
i =1 

where weights Wi can be interpreted as the total number of each stock available in the market. 

Since a comonotonic random vector X c := (Xf, ■ ■ ■ , X'f) can be represented as 

the market index under the comonotonic stock prices assumption, assuming the marginal distributions of 
individual stock prices to be unchanged, can be defined as 

d 

S c :=J2 w iFx! (V) 

1=1 

Noting the fact that, as shown in Remark 3, the variance of the market index is maximized when the 
individual stock prices are comonotonic, the herd behavior index by Dhaene et al. (2012) is defined as the 
ratio of variance of the market index to that of the index under the comonotonic assumption. The following 
definition defines the simplified version of HIX. The original version of HIX defined using the option prices 
can be found in Dhaene et al. (2012). 


HIX 



Var[ff] 
Var[5 c ] ’ 


While HIX is a convenient measure which can measure the herd behavior effectively, HIX may be sensitive 
to the marginal distributions (Choi et al., 2013). The revised version of HIX (RHIX) defined as 


RHIX xT) 


WiWjcav(Xi,Xj) 

_ 

Z WiWjCOv (xf, xf\ 
iXi V J 


is known to reduce the marginal distribution effects (Choi et al., 2013; Lee and Ahn, 2014a). The same 
measure was proposed by Dhaene et al. (2014a) from a slightly different perspective. 

Importantly, original definition of HIX (hence RHIX) can be calculated using the individual option 
prices and option price of the market index (Dhaene et al., 2012; Linders and Schoutens, 2014), and these 
measures can be used as predictors of the degree of herd behaviors in the future as implied by current option 
prices. These are the main reasons why HIX and RHIX arc favorable herd behavior indices, although there 
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may be some preference between HIX and RHIX. Of course, HIX and RHIX can also be estimated from the 
high frequency stock market data (Lee and Ahn, 2014a). 


6.2. Marginal Dependency of RHIX 

Despite some controversy, if the perfect herd behavior corresponds to comonotonic movement (Dhaene 
et al., 2012), the herd behavior should be a phenomenon that solely depends on the copula. In this sense, 
RHIX may be more favorable than HIX, because it is known to reduce the marginal distribution effects 
(Choi et al., 2013). However, as expected from the definition of RHIX (it is defined based on covariances), 
RHIX cannot thoroughly remove the marginal effects. Through a simple example, this section explains that 
such marginal effects in the calculation of RHIX can be arbitrarily large. 

For expository purposes, we consider the following Toy Model using the bivariate lognormal distribu¬ 
tion, which is frequently used to describe the stock prices. 

Toy Model. Consider only two assets 2? = {Xi.Xf) that follow a bivariate lognormal distribution with 
drift vector it and covolatility matrix X which is defined as 

s = ( ^ pa T 1 . 

V pc 1(7 2 Cf 2 ) 

Note that it. c\ and oy are parameters related with marginal distributions, and p is the only parameter for 
the (Gaussian) copula; see, for example, Nelsen (2006) and Cherubini et al. (2004) for more details. 

Under the Toy Model, simple calculation shows that 

RHIX(w,Xm eXp( ( 12J1 T l ~ 1 - (31) 

exp(cri<j 2 ) - 1 

We refer to Choi et al. (2013) for more detailed calculation for HIX and RHIX under the lognormal model. 
Now from (31), the following equality shows that RHIX under the Toy Model converges to 0, when the 
common volatility coefficient ay = ay increases, regardless of the copula coefficient p. 

lim RHIX 2^) =0 for any p < 1. (32) 

<71=172 —.>00 V / 


Knowing the degeneracy of RHIX as described in (32), the convergence rate can be important for the prac¬ 
tical use of RHIX. The following example shows the convergence rate of RHIX in (32) under various cir¬ 
cumstances. 


Example 2. Figure 2. (a) and (b) show the variation of RHIX in the Toy Model depending on the varying 
common volatility a := c\ = oy on the different scale time intervals of (0, 0.5) and (0,5) respectively. 
As shown in Figure 2. (a), RHIX looks stable around reasonable weekly volatilities of the stock markets 
assuming the weekly volatility to be 0.03. 1 RHIX around yearly volatility (= 0.03 • s/52 ~ 0.22) even 
looks stable. However, Figure 2. (b) shows that RHIX slowly but surely decreases and converges to 0 as a 
increases. 


'The weekly volatility for the S&P 500 index and IBM from March to May of 2003 were 0.0309 and 0.0365 respectively. 
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(a) RHIX for various p on the interval (0, 0.5) 


(b) RHIX for various p on the interval (0, 5) 


Figure 2: RHIX with volatility effects. 


6.3. New Herd Behavior Index: The Marginal Free Measure of Dependence 

In the following definition, we propose a new herd behavior index that is free of marginal distribution 
and hence defined in terms of copula only. 

Definition 5. For a given random vector X, Spearman’s rho type of the Herd Behavior Index (SIX) is 
defined as 

E WiWjP2 (. Xi,Xj ) 


SIX 


(tit,?) := 


i<j 


E WiWjP2 [Xf,X(j 

i<j ' 

E WiWjP2 (■ Xi,Xj ) 

i<j 


E W{Wj 

i<j 

where Spearman’s rho p 2 is defined as 

P2(X i ,X j ) = 3 (P ((Xi - X*)(Xj - X**) > 0) - P ((Xi - X*)(X, - X**) < 0)) 

with (X*,X*) and (X**,X**) are independent copies of (2Q, Xj). Sometimes we use SIX (lit, IF) to 
denote SIX (l3, X j. Note that SIX coincides with pair-wise Spearman’s rho defined in Schmid and Schmidt 
(2007) with the equal weights w\ = ■ ■ ■ = w & 

Since bivariate Spearman's rho does not depend on the marginal distribution, clearly SIX does not 
depend on the marginal distributions. Hence, for continuous marginals, we have 

SIX (Ft, FT) = SIX (Ft, C). 

Since SIX can be obtained by replacing the covariance terms in RHIX with the Spearman’s rho terms, 
it can be interpreted as the ratio of the weighted pairwise Spearman’s rho of stock prices to the weighted 
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average of Spearman’s rho of stock prices under the comonotonic assumption. Furthermore, similar to RHIX 
as in Lee and Ahn (2014a), SIX can be expressed as the weighted average of the pairwise Spearman's rhos 
as shown below 

SIX (it, = E [Z] 

where 

P(Z = p 2 (X h Xj)) = pij 

with 

WiWn 

Pi, ■= d -• 

E w kWi 

k^l 

Unlike HIX or RHIX, the calculation of SIX using the vanilla option prices may be difficult in reality 
because, whereas the calculation of HIX and RHIX requires the option prices on the individual stocks and 
the market index, the calculation of SIX requires the option prices related to every pairs of the individual 
stock prices. As an alternative, high frequency stock price data can be used for the estimation of SIX: a 
detailed method for the estimation of HIX and RHIX using high frequency stock price data can be found in 
Lee and Ahn (2014a) and a similar method can be applied to the estimation of SIX. Empirical analysis of 
herd behaviors in the stock market using stock price data and SIX can be found in Lee and Ahn (2014b). 

Remark 4. In the calculation of HIX and RHIX, we have to calculate the variance or covariance under 
the comonotonic assumptions. Hence, in the calculation of HIX and RHIX, an assumption on the marginal 
distributions is essential as shown in Lee and Aim (2014a), where lognormal distributions are assumed. 
However, for the calculation of SIX, since Spearman’s rho under comonotonic assumption is always 1 
regardless of the marginal distributions, marginal assumption is not necessaty. 

The following example present the representation of SIX in the multivariate log-normal distribution, and 
confirms that SIX is free of marginal distribution. 

Example 3. For Til E (0. oo) d and d-variate log-normal random vector A = (X\ . ■ ■ ■ ,Xf) with drift 
vector ft and covolatility matrix S, SIX can be represented as 

d 

E WiWj^ aersin (pij/2) 

SIX (w, X) = ^--- 

v 7 a 

E WiWj ( 33 ) 

d g 

= 22 Qj -aersin (pij/2) 


where p t j 



and the first inequality is from Kendall and Gibbons (1990). 


Spearman’s rho preserves the concordance ordering, and one can easily expect that SIX also preserves 
the concordance ordering. Hence it is possible to show that the maximum of SIX is achieved with the 
comonotonic copula. However, due to the absence of the minimum copula in the Frechet Space, the min¬ 
imum of SIX is not as clear. The following theorem provides some properties of SIX and determines the 
maximum and minimum of SIX. 
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Theorem 3. For given w, define S\ := w\ + • • • + Wd and S 2 '■= uo\ + • • • + w^. Then , for the given 
distribution functions H := C{F\, • • • , Ff) and H* := C*(Fi, ■ ■ ■ , Ff), the following holds: 

i. If copulas C , C* € F d satisfy C < C*, then SIX (^, H) < SIX (^, if*). 

ii. SIX satisfies 

^L_[l2l(^)-s 2 ] <SlX(ti,H) < 1, (34) 

where the definition ofl(-) can be found in Notation 2. 

Hi. The upper bound of (34) is attained if and only if H is comonotonic. 

iv. The lower bound of (34) is attained if and only if H is w*-CM, where u ? is defined in (24) of Notation 
2 . 


Proof. The proof of part i comes from the concordance property of Spearman’s rho and the fact that SIX is 
a lineal - combination of bivariate Spearman’s rho. For the proof of the remaining parts, note that 


SI X(vfr,H) = SI X(^,C) 



= C 2 - c ^ Wi + ^2 w i E [Ui] + 2 ^ WiVJjE [UiUj\ 

i =1 i =1 i<j 

= + 2 ^2 WiWjCav(Ui, Uj ) 

i<j 


where a constant c is defined as 



Now, Theorem 3 and (35) derive that 


(35) 


l(v$) < yhS 2 + 2 WiWjCav(Ui, Uj) < 

i<j 



(36) 


where the first equality in the first inequality is achieved if and only if U is w> -CM, and the second inequality 
is achieved if and only if T? is comonotonic. Now (36) and the following observation 


p 2 (X i ,X j ) = p 2 (U i ,U j ) 

= 12 Cov(Ui,Uj ), 
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conclude the following inequalities 


^ [12 l(ti) - S 2 ] < SIX (1&, H) < 1 

where the first equality holds if and only if H is re*-CM and the second inequality holds if and only if H is 
comono tonic. 

□ 


6.4. Data Analysis 


In this subsection, we analyze the herd behaviors in the US stock market using SIX. Daily stock prices 
(t) of three stocks arc collected from Apple, Hewlett-Packard Company and New York Times in the 
time interval between t =2001/March/01 and t =2014/April/09. Under the lognormal model, the line 
graph (SIX) in Figure 3 shows estimated SIX, where SIX at each point is estimated based on 4 month 
observations. Similar to Choi et al. (2013), three stock prices shows generally strong herd behavior during 
the global financial crisis starting from 2008. 

Sometimes, we may be interested in the relationship between the the stock prices of three companies 
only. For example, we may assume that the stock prices of three companies reflect the preferences between 
the traditional media system (newspapers), the traditional internet based media system (computers), and 
the mobile internet based media system (smartphones or tablets). However, strong dependency of the three 
stock prices may not stand for the strong dependency between three companies in particular, because the 
strong comovement of the stock prices during the period may be the result of the illusion effect caused by 
devaluation of the whole stock market (the global economic crisis in 2008, for example). Hence, to under¬ 
stand the actual physical relation between three stock prices, it can be beneficial to consider the detranded 
stock price by the market index (S&P in this data analysis) defined as follows: 


X M (t ) := 3(t)/S(t), 


where S(t) is S&P index. The dashed line graph (SIX M )in Figure 3 shows estimated SIX using the adjusted 

stock price X AI on the same time interval. Here, we have used the weight t? = (1,1,1). Note that, under 
the lognormal model in (3), specific statistical estimation procedures can be found in Lee and Ahn (2014a), 
for example. 

From Figure 3, we can conclude that main source of the comovement during the global financial crisis 
is the devaluation of the whole stock market. After removing the comovement effect by the global financial 
crisis, comovement of adjusted stock prices 2 ? M is not as strong. 


7. Conclusion 


In this paper, we have provided the set of copulas called i^-CM copulas, and have shown these to be 
the minimal in set concordance ordering. Given the absence of the minimum copula, the minimality can 
be important in optimization problems. Especially, we show that the proposed set of copulas minimize the 
variance of the aggregated sum where the marginal distributions arc given as various uniform distributions. 
As shown in Remark 3 in Section 5, the set of minimal copulas can be related with the variance of aggre¬ 
gated sum with given marginal distributions. In this respect, the approach using d-CM copulas, which arc 
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Figure 3: SIX from 2001/March/01 to 2014/April/09 with weight vector (w\. w- 2 , v'n) = (1,1,1). 


the generalized version of 77’-CM copulas, can be shown to be useful in minimizing the variance of the 
aggregated sum for some special marginal distributions. We leave this topic for future research. 

Finally, although 7?-CM copulas do not minimize the variance of the aggregated sum in general when 
the marginal distributions are not uniform distributions, many other interesting optimization problems have 
uniform marginals as their solutions. Optimization of the herd behavior index is one such example. In 
this paper, we have provided a herd behavior index that does not depend on the marginal distributions, and 
showed that the herd behavior index is minimized with Tit-CM copulas. 
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Appendix A. 

Lee and Ahn (2014b) proposed the class of minimal copulas which can be viewed as alternatives to 

countermonotonicity in multivariate dimensions. 
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Definition 6 (Lee and Ahn (2014b)). A d-variate random vector U will be called d-countermonotonic 
(d-CM) or d-CM if there exist function (/i, • • • , fd) £ A4({_[0,1] and 

d 

Y.f'^') = c (A.1) 

3 =1 

with probability 1 for some constant c£l Equivalently, we say that the distribution function C is d-CM 
if iJ is d-CM. Especially, for the choice of functions with c = 1 in (A. l),Tf is called d-CM with parameter 
functions (/i, • • ■ , fd). 

Since Lee and Ahn (2014b) have shown that d-CM does not depend on marginal distributions (see 
Lemma 1 in Lee and Ahn (2014b)), we provide a version of d-CM definition for a copula only in this 
appendix. As we have briefly mentioned in Section 3, d-CM may be too general to be used for the extreme 
negative dependence as it includes almost countermonotonic movement. Alternatively, Lee and Ahn (2014b) 
provide a definition of strict d-CM as a subset of d-CM in the following sense. 

Definition 7 (Lee and Ahn (2014b)). A d-variate random vector iJ is strict d-CM if 


P 




= 1. 


Equivalently, we say that H is strict 



is strict d-CM. 


It is obvious that strictly d-CM is d-CM having constant multiplication of identity functions as parameter 
functions: i.e. 

fi(v) = ■■■ = f d (v) = 2' v 

for v £ [0,1]. The existence of a strict d-CM copula is shown in Ruschendorf and Uckelmann (2002); Lee 
and Ahn (2014b). Strict d-CM is useful in various minimization/maximization problems (Lee and Ahn, 
2014b). 


Proof of Corollary 1 . Showing Corollary 1 is equivalent to show that for any given 7?-CM copula C and 
C* E Fd satisfying 

C* PC, 

implies that C* is also rt^-CM. 

First observe that if Tit does not satisfy (16), then C is empty and the proof is trivial. So we can assume 
that Tit satisfies (16) and C is not empty. Now, define two sets 


F r := { 


it £ [0, l] c 


d 

E w i 

1=1 


d 

WiUi < 2 

i= 1 
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and 


Qc '■= < 


it e [o,i] f 


d 

d Y, w i 

J2wm < —^— and ?ii, 

i =1 


, Ud arc rational numbers 


Then, by the denseness of rational numbers in real line, we have 


it G [0, l] c 


d 

d J2 w i 

Y w i x i < 

i= 1 


U {^e[0.l]“ 

1 t£Fc 

U {^e[0,i]‘ 

i teQc ^ 


it < it 


it < it > . 


which in turn implies 

( 

P 


E w i 

Ywi.U*<^— 


i —1 


v 


\ 

/ 

/ 

= p 

U * G < 

it € [0, l] d 

/ 

V 

\ 


* )\ 

E w 


Y < i ~y 


i= 1 


p (tr e [J g [o,i] d E < 


TteQ, 

< y p(p <^) 

iteQc 

= 0. 

where the last inequality holds because Q c is countable set. Si mi lar logic derives 

/ 


Y w * u ? 


> 


d \ 

E wi 

1 


i= 1 


v 


= 0 


which in turn concludes the proof with (A.2). 


(A.2) 


□ 
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